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EULERIAN  CURRENT  MEASUREMENTS  AT  PHELPS  BANK 


INTRODUCTION 

From  July  5-25,  1982  scientists  from  NRL  and  other  laboratories  con¬ 
ducted  a  cooperative,  remote  sensing  experiment.  The  general  purpose  of 
the  experiment  was  to  obtain  information  on  oceanographic  processes 
responsible  for  the  surface  expression  of  bathymetry  (SEBEX)  and  hydro¬ 
graphy  in  the  wave  field  and  radar  imagery  of  relatively  shallow  seas.  To 
this  end,  simultaneous  and  coordinated  remote  sensing,  oceanographic, 
meteorological,  hydrographic  and  bathymetric  measurements  were  made.  The 
site  chosen  as  the  central  point  for  the  experiment  was  Phelps  Bank.  The 
bank  is  a  relatively  isolated,  subsurface,  topographic  feature  located 
approximately  37  nautical  miles  southeast  of  Nantucket  Island 
( 40°50 1 N-69°20 fW) .  Figure  1  shows  where  the  bank  appears  on  the  navi¬ 
gational  chart  (NOAA,  1979). 

The  field  exercise  was  under  the  direction  of  Dr.  Davidson  Chen  (NRL 
Code  7912C)  and  Dr.  Gaspar  Valenzuela  (NRL  Code  4305).  Mr.  William 
Garrett  (NRL  Code  4333)  served  as  senior  scientist  aboard  the  USNS  HAYES. 
The  overall  plan  of  the  NRL  Remote  Sensing  Experiment  has  been  submitted 
for  publication  by  Valenzuela  (1981)  and  a  review  of  surface  effects 
attributable  to  subsurface  processes  has  been  published  by  Chen  (1982). 
This  report  summarizes  the  Eulerian  current  measurements  made  during  the 
remote  sensing  experiment. 

BACKGROUND 

Among  the  first  investigators  to  report  that  radar  imagery  of  the 
surface  waves  of  shallow  seas  included  features  that  were  related  to  the 
Manuscript  approved  January  19,  1983. 
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bottom  topography  were  DeLoor  and  Brunsveld  van  Hulten  (1978),  The 
phenomenon  has  been  treated  in  some  detail  in  DeLoor* s  more  recent  work 
(DeLoor  1981).  The  bottom  topography  effect  has  been  particularly  notice¬ 
able  in  SAR  (synthetic  aperture  radar)  imagery  from  the  SEASAT  satellite 
(English  Channel  and  Nantucket  Shoals)  although  airborne  SAR  and  SLAR 
(side-looking  airborne  radar)  also  show  the  effect.  McLeish  et  al  (1981) 
have  published  an  interpretation  of  SLAR  imagery  of  surface  wave  patterns 
in  the  North  Sea  which  are  related  to  sea-floor  topography.  General 
reviews  of  radar  imagery  of  the  ocean  surface  may  be  found  in  the  recent 
book  "Spaceborne  synthetic  aperature  radar  for  oceanography"  Beal  et  al. 
eds.  (1981)  and  in  Alpers  et  al.  (1981). 

In  order  to  interpret  surface  manifestations  of  subsurface  topography 
it  is  necessary  to  have  information  regarding  the  flow  of  current  over  and 
around  the  particular  topographic  feature  (in  this  case  Phelps  Bank). 
Currents  at  the  site  were  anticipated  to  be  predominantly  rotary  tides  (1 
to  2  knots),  based  on  data  provided  for  the  nearest  locations  in  the  NOAA 
Tidal  Current  Tables  published  by  the  U.S.  Department  of  Commerce  (NOAA, 
1981).  The  locations  provided  by  the  tide  tables  are  Nantucket  Shoals 
40°37tN-69°37,W;  Davis  Bank  East  41°02 f N-69°41 rW;  and  Great  South  Channel 
40°31fN-63°47  fW=,  These  locations  are  13,  19  and  32  nautical  miles 
respectively  fr;m  Phelps  Bank.  The  tidal  predictions  at  the  three  sites 
are  based  on  the  time  of  maximum  flood  tide  at  Pollock  Rip  which  is 
located  approximately  48  nautical  miles  northwest  of  Phelps  Bank.  Pro¬ 
gressive  vector  diagrams  of  the  current  speeds  and  directions  were  plotted 
for  the  three  sites.  The  major  axes  of  the  resultant  tidal  ellipses 
varied  in  direction  from  345°  to  40°  and  the  current  speeds  ranged  between 
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0.8  and  1.7  kt.  for  the  same  tidal  phase.  This  relatively  large  variation 
in  both  current  speed  and  direction  at  the  locations  provided  by  standard 
tide  tables  clearly  indicated  the  need  for  real-time  current  measurements 
at  the  specific  site  in  question  (Phelps  Bank). 

FIELD  MEASUREMENTS 

Currents  in  the  vicinity  of  Phelps  Bank  were  measured  by  an  Eulerian 
method  in  addition  to  the  Lagrangian  drogue  measurements  that  have  been 
reported  previously  (Greenewalt  and  Gordon,  1982).  The  Eulerian  tech¬ 
nique  consists  of  recording  the  speeds  and  directions  of  currents  as  they 
flow  past  a  fixed  point.  In  this  case  the  fixed  point  was  a  moored 
array  of  three  recording  current  meters.  The  mooring  was  located  at  a 
position  approximately  1.2  naut.  mi.  east  of  the  shallowest  part  of  Phelps 
Bank  (40°  50.07 fN  -  69°  l9.81fW  as  marked  in  Figure  1.).  The  depth  at  the 
site  is  nominally  30  m  with  an  uncertainty  estimated  to  be  ± 2  m.  The 
uncertainty  accounts  for  errors  in  depth  measurement  and  tide  level  range. 
The  mooring  system  is  shown  schematically  in  Figures  2  and  3.  The 
buoyancy  (Fg  approximately  110  lb.)  was  provided  by  two  1?”  diameter 
hollow  glass  spheres.  This  flotation  supported  a  1/4"  diameter  steel 
cable  which  was  attached  to  a  250  lb.  anchor.  A  surface  buoy  to  aid  in 
relocation  and  recovery  of  the  mooring  was  connected  to  the  subsurface 
flotation  by  150 f  of  polypropylene  line.  The  recording  current  meters 
were  attached  to  the  steel  cable  at  the  positions  shown  in  Figure  2. 

Under  conditions  of  slow  currents  the  meters  designated  at  M^,  M2,  and  M3 
were  at  nominal  depths  of  20.6,  11.4  and  3.4  meters  respectively.  In  the 
presence  of  the  stronger  currents  the  actual  depths  of  the  current  meters 
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Fig.  2.  -  A  schematic  diagram  of  the  forces  acting  on  the  current  meter 


mooring. 


CURRENT  METER  MOORING 


to  the  bathymetry  of  Phelps  Bank.  The  vertical  scale  is 
exaggerated  by  a  factor  of  X36. 


Table  1.  The  effect  of  drag  on  the  mooring 


CURRENT 

<KT) 

DRAG 

(Fd  LBS.) 

DISPLACE¬ 
MENT  ANGLE 
(0  DEG) 

|  CURRENT  METER  DEPTH  (M) 

mm 

m2 

_ 

m3 

0.25 

1.1  j 

0.5 

- j 

i 

20.63 

mmgm 

3.41 

0.50 

4.4  1 

2.2 

20.64 

3.43 

0.75 

10.0 

5.1 

20.67 

3.51 

1.00 

17.7 

9.0 

20.74 

3.73 

1.25 

27.7 

13.9 

20.90 

1 

4.19 

1.50 

39.9  : 

19.6 

21.17 

'SIB 

4.95 

1.75 

54.3 

25.9 

21.57 

1 

6.07 

2.00 

70.9 

32.3 

22.08 

14.29 

7.52 

2.25 

89.8 

1 _ i 

38.7 

22.69 

15.50 

9.26 

It  should  be  noted  that  in  Table  1  the  significant  figures  for  meter  >ths 

represent  only  the  geometry  and  not  absolute  precision*  It  is  seen  f 

the  table  that  under  the  given  conditions  meter  M3  may  vary  in  depth 

as  much  as  4  meters  in  the  current  range  from  one  to  two  knots*  Below  one 

knot  the  meters  change  depth  very  little.  This  is  a  consequence  of  the 

fact  that  the  drag  force  is  proportional  to  the  square  of  the  current, 

9 

i.e.,  Fry-V .  The  major  sources  of  uncertainty  in  the  calculations 
for  Table  1  are  the  drag  coefficient  Cq  and  the  assumption  that  the  drag 
is  localized  at  the  flotation.  The  overall  error  is  estimated  to  be  the 
order  of  ±20%. 

The  current  meters  attached  to  the  mooring  are  of  the  pressure-vane 
type,  that  is,  their  suspension  angle  depends  on  the  current  roughly  in 
proportion  to  .  The  devices  are  film  recording  current  meters,  manu¬ 
factured  by  General  Oceanics,  Inc.  of  Miami,  Florida.  Figure  4  is  a  cut¬ 
away  drawing  of  the  meter.  In  the  words  of  the  manufacturer  it . 

'* . consists  of  a  finned  cylindrical  housing  containing  a  direc¬ 

tional  inclinometer  and  Super-8  cartridge  camera  which  sense  and 
record  the  inclination  and  compass  heading  of  the  instrument. 
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A.  -  A  cut-away  drawing  of  the  configuration  of  General  Oceanics , 


vane  current  meter# 


Fins  are  affixed  to  the  current  meter  housing  to  assist  in  direc¬ 
tional  orientation  and  stabilization.  The  fin-housing  combination 
within  the  current  stream  creates  a  drag  resulting  in  an  inclina¬ 
tion  of  the  instrument  from  the  vertical  which  changes  in  direct 
relation  to  the  current  velocity. 

The  data  recording  camera  is  triggered  to  photograph  the 
inclinometer  by  a  quartz  crystal  controlled  timing  circuit  con¬ 
tained  within  the  camera.  The  self-contained  battery  supply  and 
the  film  cartridge  capacity  enable  up  to  7200  data  records  to  be 
taken  over  operating  periods  of  up  to  ten  months. 

The  directional  inclinometer  is  a  spherically  shaped  compon¬ 
ent  mounted  on  the  inner  face  of  the  lower  housing  end  cap.  The 
inclinometer  design  utilizes  a  transparent  fluid  filled  housing 
containing  a  negatively  buoyant  inner  sphere  floating  on  a  bear¬ 
ing  of  liquid  mercury.  The  inner  sphere  maintains  a  stable  vert¬ 
ical  attitude  and  magnetic  north  orientation  because  of  an 
internal  bar  magnet  whose  mounting  location  gives  the  sphere  a  low 
center  of  gravity.  A  small  circular  target  at  the  top  of  the 
transparent  housing  is  viewed  by  the  camera  against  a  grid  of 
precision  latitude  and  longitude  lines  inscribed  on  the  free 
inner  sphere.  When  photographed  by  the  camera,  this  target  mark 
enables  direct  reading  of  the  instrument  attitude  and  azimuth  by 
its  position  relative  to  the  latitude  and  longitude  lines. 

Also  located  on  the  lower  end  cap  above  and  to  one  side  of 
the  directional  inclinometer  is  a  battery  powered  quartz  crystal 
watch  which  provides  the  time,  day  and  date  for  each  data  frame." 

heedless  to  say  the  data  interpretation  procedure  is  rather  tedious, 

that  is,  each  frame  of  motion  picture  film  provides  only  one  reading  of 

current  speed  and  direction.  The  directional  inclinometer  is  read  as 

follows : 


"The  inner  sphere  of  the  directional  inclinometer  is  gradu¬ 
ated  in  increments  of  10°.  Every  circle  of  latitude  is  equal  to 
10°,  starting  at  zero  and  finishing  at  90°.  The  stamped  circles, 
i.e.,  3,  3,  7,  &  9,  indicate  30°,  50°,  70°  and  90°  respectively 
from  the  vertical. 

The  lines  of  longitude  are  graduated  in  thirty-six  10° 
increments  for  a  total  of  360  degrees.  Every  letter  indicates 
30°  of  compass  heading,  e.g.,  ,#D"  at  060°;  MN"  at  0°  (or  360°)." 

Figure  5  shows  how  the  circular  target  appears  between  the  lines  of 

inclination  and  magnetic  orientation.  In  Figure  5  the  deflection  is  65° 
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and  the  magnetic  orientation  is  225°,  i.e.,  between  south  and  west.  The 
65°  tilt  or  deflection  angle  is  related  to  the  current  speed  and  can  be 
converted  to  speed  using  the  calibration  curves  provided  by  the  manufac¬ 
turer.  Figure  6  is  an  example  of  such  a  calibration  curve  and  using  it 
the  65°  deflection  is  found  to  be  equivalent  to  a  current  speed  of  52 
cm/sec.  It  should  be  noted  that  during  the  measurements  reported  here 
the  currents  were  relatively  fast  (1-2  knots)  and  the  deflection  angles  of 
the  current  meters  were  near  the  maximum  of  their  calibrated  range. 

Because  of  this  and  variations  in  the  mountings,  the  absolute  accuracy  of 
the  measurements  should  be  considered  no  better  than  ±10%. 

For  the  measurements  described  here  the  current-meter  mooring  was 
deployed  at  about  0910  local  time  July  10,  1982  and  recovered  at  1800 
local  time  on  July  20.  The  duration  of  the  measurement  time  series  was, 
therefore,  249  hours.  The  cameras  recorded  8  frames  per  hour  in  the 
three  current  meters  for  a  total  of  about  6000  readings  of  current  speeds 
and  directions.  All  6000  readings  are  not  listed  here  in  tabular  form, 
however  all  the  values  have  been  stored  on  tape  and  can  be  made  available. 
The  following  tables  (Tables  2-4)  list  hour  by  hour  averages  of  the  current 
speeds  and  directions  (8  data  points)  for  the  10  days  the  mooring  was 
deployed.  The  current  averages  include  one  half  hour  before  and  one  half 
hour  after  the  times  given  in  the  tables.  For  purposes  of  comparison  with 
the  current  meter  arrangement  shown  in  Fig  2.  Meter  "X"  was  at  a  nominal 
depth  of  about  5  m  and  corresponds  to  M3;  "M"  corresponds  to  M2  at 
about  13  m  and  "NM  corresponds  to  M]^  at  21  m.  The  times  in  the  tables 
are  in  universal  time  (U.T.).  To  convert  to  local  time  (EST)  subtract 
four  hours  or  to  obtain  Daylight  Saving  Time  (EDST)  subtract  5  hours. 
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ion  curve  for  a  General  Oceanics  vane  current  meter. 
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Table  2.  Current  Measurements  at  Phelps  Bank  (40°  50.07’^ 


69°  19.81fVJ) 


TIDE  Tft&LE  DERIVED  FROM 
CURRENT  METER  1 X ‘  ft T  DEPTH  OF 
FIVE  METERS 

Speeds  and  directions  are 
averages  o t  e i e h t  measure m e n t 5 
The  table  aives  hourly  speeds 
and  directions  tor  the  ten  day 
period  from  July  10  to  20.*  1982 
Hours  are  in  universal  time. 
Subtract  tour  hours  tor  local 
time  or  five  hours  for  Eastern 
S  t  andard  Time. 


DRY  l SI  JULY  10  1582 

HOUR  SPEED  DIRECTION 
u t  k t s .  decrees 


15 
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X 

98 
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5 

48 
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1 

48 
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1 
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11  193 
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1  46 
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1  .  80 

L  ~  i  t7 

77 

.  64 

30  5 

4 
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5 
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6 
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7 
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26 
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1 38 
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2  2 
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HOUR 

SPEED 

DIRECT 

u  t 

kts . 

decree 

0 

1 .41 

236 

1 

1 . 30 

246 

2 

1  .  17 

257 

7 

1.01 

267 

4 

.  53 

333 

5 

.  54 

9 

K 

.  54 

70 

( 

.  74 

126 

Q 

.  90 

1  70 

Q 

1  .  27 

199 

10 

1 . 38 

219 

1  1 

1  .  65 

2  3  y 

12 

1 . 52 

237 

13 

1 . 36 

248 

14 

1.19 

2  72 

15 

1.07 

386 

16 

85 

356 

17 

1  .  15 

45 

18 

1  .  40 

(  i 

19 

1 . 38 

95 

20 

1  .  23 

121 

21 

1.31 
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2  2 

1  .  32 
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2 

1 . 32 
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i 
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cr  7 
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cr  o 
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6 

67 

79 

i 

.  88 
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x 
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177 

9 

1 
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10 

1 

!  65 
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1 
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12 

1 

65 
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I  3 

1 

79 
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14 

1 
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15 

75 

3 1 8 

16 

1 

02 

C1 
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1 
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54 

18 

1 
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l‘  w 

19 

1 

er 
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20 

1 
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21 

1 
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1 
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Table  3*  Current  Measurements  at  Phelps  Bank  (40°  50.07 fN  -  693  19.81’W) 


TIDE  TABLE  DERIVED  FROM 
CURRENT  METER  1 M '  AT  DEPTH  OF 
THIRTEEN  METERS 
Speeds  and  directions  are 
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Table  4.  Current  Measurements  at  Phelps  Bank  (40°  50.07'N  -  69’  19.81'W). 
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Table  4.  (Continued) 
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Two  qualitative  aspects  of  the  currents  can  be  established  from  simple 
inspection  of  the  tables.  First  of  all,  it  is  seen  that  the  current  direc¬ 
tion  varies  continuously  in  a  clockwise  manner  at  a  rate  of  approximately 
29°  per  hour.  That  is,  the  current  at  Phelps  Bank  is  dominated  by  a  rotary 
tidal  component  with  a  period  of  about  12-1/2  hours.  It  is  also  seen  from 
the  tables  that  the  tidal  current  is  fairly  strong,  occasionally  reaching 
values  of  1.8-1. 9  knots  (.92-. 98  msec"-*-).  The  slowest  currents  fall  in 
the  range  of  0.50  knots  (.26  msec“^).  The  high  currents  are  associated 
with  directions  in  the  210°-260°  range  while  the  slow  currents  are  corre¬ 
lated  with  directions  between  290°  and  about  110°.  This  indicates  that 
the  trajectory  of  the  water  during  a  tidal  cycle  is  generally  elliptical. 

The  information  in  Tables  2-4  is  presented  graphically  in  Fig.  7,  A-C , 
which  shows  the  hourly  averaged  values  of  current  speed  and  direction  for 
all  three  meters  during  the  10  days  the  mooring  was  in  place.  Perhaps 
the  most  striking  feature  of  Fig.  7  is  that  the  current  meters  at  diff¬ 
erent  depths  show  no  gross  differences  in  their  speeds  and  directions. 

For  example,  the  speed  maxima  and  minima  have  very  nearly  the  same  values. 
The  most  noticeable,  consistent  variation  is  that  the  direction  of  the 
deep  meter  (N)  almost  always  lags  behind  the  upper  two  meters  by 
approximately  30°.  The  only  other  obvious  consistency  is  that  the  current 
speed  measured  by  the  deep  current  meter  is  always  about  1/4  kt.  slower 
than  the  upper  meters  when  the  flow  direction  is  in  the  100°-140°  range. 

The  physical  mechanisms  responsible  for  these  effects  are  not  clear  but 
they  may  be  related  to  veering  due  to  bottom  drag  or  to  local  topographic 
inf luence . 
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urrent  meters,  July  10-20,  1982  at  Phelps  Bank. 


Current  measurements  from  a  fixed  mooring  can  also  be  represented 
graphically  as  progressive  vector  diagrams.  This  is  the  traditional 
method  for  producing  a  current  "quasi-trajectory"  based  on  measurements 
made  at  a  fixed  site.  This  is  accomplished  by  considering  each  measured 
speed  and  direction  as  a  vector  and  sequentially  adding  these  vectors  to 
each  other  in  a  time  series.  An  example  of  this  approach  is  shown  in 
Fig.  8.  Here  all  the  current  plus  direction  records  (8  per  hour)  from 
the  near-surface  current  meter  X  (~  5  m  depth)  are  vectorially  added  for 
the  10  day  duration  of  the  deployment. 

When  the  data  are  presented  in  this  way,  two  additional  aspects  of 
the  current  at  this  location  become  evident.  The  first  and  most  obvious 
is  that  there  is  a  progressive  drift  toward  the  southwest  (220°),  that  is, 
the  current  includes  a  relatively  steady  component  in  addition  to  the 
rotary  tidal  motion  mentioned  earlier.  The  second  noticeable  feature  is 
that  the  relative  amplitudes  of  sequential,  12.5-hour  rotary  tides  change 
along  the  time  series  from  July  10  to  July  20  (Julian  day  191  to  201). 

This  pattern  is  typical  of  semi-diurnal  tides  and  is  associated  with  the 
tilt  of  the  earth's  axis  of  rotation  relative  to  the  ecliptic,  the  lati¬ 
tude  of  the  observations  and  the  time  in  the  lunar  month.  Detailed  dis¬ 
cussions  of  this  tidal  phenomenon  are  available  in  most  standard  oceano¬ 
graphic  reference  texts  (Hansen,  1962)  and  it  will  not  be  pursued  here. 

We  will  now  return  to  a  consideration  of  the  long-term  flow  over 
Phelps  Bank,  It  is  seen  from  Fig,  8  that  there  is  some  variability  in  the 
speed  and  direction  of  the  drift  current  with  an  average  speed  of  approxi¬ 
mately  0.47  kt  (22  cm  sec“^).  It  should  be  noted  here  that  the  other 
two  current  meters  at  13  (M,  M2)  and  21  m  (N,  M]^)  depth  do  not  record  the 
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same  "steady"  component  as  meter  X  (M3)  at  5  m,  which  is  shown  in  Fig.  3. 
The  relative,  quasi-trajectories  of  the  steady  components  of  current  as 
recorded  by  the  3  meters  are  illustrated  in  Fig.  9  and  listed  in  Table  5. 
It  can  be  seen  that  the  drift  current  indicated  by  the  mid-depth  meter 
(M,  13  m)  is  in  approximately  the  same  direction  (220°)  as  that  recorded 
by  the  near-surface  meter  (X,  5m)  but  it  is  slower  by  about  .12  kt.  The 
deep  meter  (N,  21  m)  records  the  same  current  speed  as  the  mid-depth 
meter  but  the  direction  of  the  flow  differs  by  about  10°  (210°).  Whether 
this  latter  effect  is  due  to  veering  as  a  result  of  bottom  drag,  local 
topography  or  instrumental  error  is  not  demonstrable  from  the  limited  data 
available.  What  is  clear  is  that  there  is  a  considerable  vertical  shear 
in  the  steady  drift  current.  The  shear  between  meters  X  (5  m)  and  M  (13 
m)  is  approximately  8.5  X  10“* 3  Sec~^  and  is  primarily  a  result  of  speed 
differences.  Between  meters  M  (13  n)  and  N  (21  m)  the  shear  is  about 
3.2  X  10*^  sec“^  and  is  mainly  attributable  to  differences  in  direction. 
These  shears  are  relatively  large  but  not  atypical  of  those  measured  in 
the  upper  ocean. 
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Table  5,  Drift  current  speeds  and  directions  recorded  by  the  3  current 
meters  for  the  25  hour  intervals  shown  in  Figure  9 


The  daily,  averaged  drift  currents  listed  in  Table  5  can  be  sub¬ 
tracted  from  the  total  currents  to  obtain  the  rotary  tidal  component. 

The  result  for  the  near-surface  current  meter  (X,  5  m)  is  shown  in  Fig, 

10,  The  plot  is  disconcertingly  reminiscent  of  the  classical  "NIDUS  RATTI" 
but  does  serve  to  illustrate  the  rotary  tidal  flow  and  makes  it  possible 
to  produce  tide  tables  for  the  specific  location.  This  is  accomplished  by 
averaging  the  speeds  and  directions  of  the  tidal  currents  and  expressing 
them  in  terms  of  the  time  after  maximum  flood  tide  at  Pollock  Rip.  The 
tidal  current  at  the  Phelps  Bank  mooring  site  can  then  be  predicted  for 
any  time  in  the  future  by  using  the  Pollock  Rip  data  from  published  tide 
tables  (NOAA,  1981  or  equivalent)  and  Table  6  given  here.  The  procedure 
for  obtaining  past  or  future  tides  at  various  locations  on  Nantucket  Shoals 
are  described  in  detail  in  the  National  Ocean  Survey,  Tidal  Current  Tables 
(NOAA,  1981). 
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( *>  m)  during  the  deployment.  This  represents  the  progressive 
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Kig.  14.  -  The  rotary  tidal  component  oi  the  curren'  at  meter  N  (21  m) 


III.  Meter  N,  nominal  depth  21  m. 
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The  progressive  averaging  of  the  data  in  Fig.  10  over  25-hour  and 
12.5-hour  tidal  cycles  is  shown  graphically  in  Figs.  11  and  12  respect¬ 
ively.  It  is  seen  that  the  12.5  hour  averaging  required  produce 
Table  6  removes  the  sequential  asymmetry  associated  with  the  relative 
earth-moon  positions  that  was  mentioned  earlier.  The  equivalent  12.5- 
hour  averaged  rotary  tides  for  the  deeper  meters  (M  at  13  m,  N  at  21  m) 
are  plotted  in  Figs.  13  and  14.  A  comparison  of  Figs.  12,  13  and  14 
indicates  that  the  rotary  tidal  component  of  the  current  is  almost 
identical  at  depths  of  5  and  13  m  while  the  east-west  tidal  amplitude  is 
somewhat  reduced  at  21  m  depth.  This  difference  shows  that  there  is  some 
shear  in  the  rotary  tidal  component  of  the  current  as  well  as  in  the 
"steady”  component  discussed  earlier.  At  13  m  and  shallower  the  rotary 
component  moves  essentially  as  a  slab.  Below  13  m  the  east-west  current 
is  somewhat  slower.  The  data  set  is  not  adequate  to  attribute  this  effect 
to  any  specific  hydrographic  or  hydrodynamic  cause. 
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COMPARISON  TO  LAGRANGIAN  METHOD 


When  reviewing  and  interpreting  the  information  presented  in  this 
report,  it  is  instructive  to  compare  and  contrast  these  results  with  our 
previously  published  results  of  current  measurements  in  the  vicinity  made 
by  Lagrangian  methods  (Greenewalt  and  Gordon,  1982),  The  two  reports 
should  be  considered  complementary*  The  most  striking  distinction  between 
currents  measured  on  the  bank  by  Eulerian  methods  and  those  measured  east 
and  west  of  the  bank  by  drogue  following  (Lagrangian  technique)  is  graph¬ 
ically  illustrated  in  Fig*  15*  The  plot  shows  the  trajectory  of  a  current¬ 
following  drogue  (A)  at  6  m  depth  at  a  location  centered  about  4.5  nautical 
miles  west  of  Phelps  Bank  compared  to  the  "quasi-trajectory"  (progressive 
vector  diagram)  of  currents  measured  at  the  same  time  by  Eulerian  technique 
at  the  mooring  (Meter  X).  The  meter  array  was  deployed  near  the  east  side 
of  the  bank  (nominal  depth  =  5  m,  July  10-11,  1982),  Both  tracks  show  the 
elliptical  shape  typical  of  rotary  tides,  however,  the  ellipse  measured 
off  the  bank  where  the  sea  was  about  40  m  deep  has  a  major  axis  in  a 
roughly  North-South  direction  while  for  the  ellipse  measured  over  the  bank 
(~  30  m  depth  or  less)  the  East-West  dimension  is  extended.  The  result 
shown  in  Fig,  15  is  a  typical  one  and  applies  to  all  cases  when  there  were 
simultaneous  current  measurements  by  Eulerian  methods  at  the  bank  and 
Lagrangian  methods  in  the  vicinity  of  the  bank.  We  have  interpreted  this 
enhancement  of  current  flow  across  the  short  dimension  of  Phelps  Bank  as  a 
topographic  effect  of  the  bank  itself,  namely,  an  acceleration  of  the  flow 
due  to  the  vertical  reduction  of  its  "channel  depth".  This  interpretation 
is  also  consistent  with  current  speed  changes  across  the  hank  measured  by 
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using  the  ship  (USNS  HAYES)  as  a  Lagrangian  drifter.  Whether  the  inter- 
pretation  can  be  quantitatively  confirmed  depends  on  more  detailed 
measurements  of  the  local  bathymetry  and  current  speeds  across  the  topo¬ 
graphic  feature.  Elementary  continuity  analysis  of  the  current  speed  and 
depth  measurements  across  the  bank  (Greenewalt  and  Gordon,  1982)  indicate 
that  the  situation  may  be  somewhat  more  complex  than  simple  channel  con¬ 
straint,  i.e.,  the  product  of  depth  and  current  speed  is  not  constant  over 
Phelps  Bank. 

SUMMARY 

The  primary  purpose  of  this  report  was  to  supply  the  participants  in 
the  NRL  Remote  Sensing  Experiment  (July  5-25,  1982)  with  hour  by  hour 
measurements  of  current  at  Phelps  Bank  to  aid  in  their  interpretation  of 
hydrographic,  radar,  wave,  IR,  photographic  and  other  simultaneous 
measurements  in  the  vicinity.  A  second  objective  was  to  provide  a  means 
of  predicting  currents  at  Phelps  Bank  at  some  future  time  for  purposes  of 
operational  planning.  Both  these  tasks  have  been  accomplished  with  the 
data  listings  and  tide  tables  produced. 

Comparison  to  current  measurements  made  by  Lagrangian  methods  indi¬ 
cates  that  the  "channel  depth  constraint"  of  Phelps  Bank  causes  the  East- 
West  component  of  the  rotary  tidal  currents  to  be  enhanced  as  it  flows  over 
the  topographic  feature. 

For  purposes  of  planning  future  field  exercises  of  a  similar  nature, 
it  is  recommended  that  current  meters  be  moored  on  both  sides  of  Phelps 
Bank  as  well  as  on  the  central  crest  of  the  bank.  Current  meters  should  be 
deployed  in  at  least  6  locations:  2  east,  2  west  and  2  on  top  of  Phelps 
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Bank.  The  measurements  by  Eulerian  methods  must  be  accompanied  by  simul¬ 
taneous  Lagrangian  measurements  for  purposes  of  mutual  calibration  and 
interpretation  of  the  hydrodynamics  of  the  flow  regime. 

ACKNOWLEDGMENT 

The  authors  express  their  appreciation  to  Mr.  John  H.  Ostrander,  NRL 
navigator,  and  the  NRL  Ship  Facilities  Group  for  their  invaluable 
assistance  during  deployment  and  recovery  of  the  mooring. 


42 


JJJLtLUJLSBt 


REFERENCES 

Alpers,  W.R.,  D.3.  Ross  and  C.L.  Rufenach  (1981).  On  the  Detectability  of 
Ocean  Surface  Waves  by  Real  and  Synthetic  Aperture  Radar,  J.  Geophys.  Res. 
36,  No.  C7,  6481-6498. 

Beal,  R.C.,  P.S.  De  Leonibus  and  I.  Katz,  editors  (1981).  Spaceborne 
synthetic  aperture  radar  for  oceanography,  Johns  Hopkins  University  Press, 
Baltimore. 

Chen,  D.T.  (1982).  Surface  Effects  Due  to  Subsurf aace  Processes:  a 
Survey,  NRL  Memorandum  Report  4727,  January  15,  1-40. 

De  Loor,  G.P.,  and  H.W.  Brunsveld  van  Hulten  (1978).  Microwave  Measure¬ 
ments  over  the  North  Sea,  Boundary^Layer  Meteorology,  13 ,  119-131. 

De  Loor,  G.P.  (1981).  The  Observation  of  Tidal  Patterns,  Currents  and 
Bathymetry  with  SLAR  Imagery  of  the  Sea,  IEEE  Jour,  of  Oceanic  Eng. 

0E-6,  No.  4,  124-129. 

Greenewalt,  D.  and  C.  Gordon  (1982).  Lagrangian  Current  Measurements  at 
Phelps  Bank,  NHL  Memorandum  Report  4965,  57  pp,  December  6,  1-55. 

Hansen,  W.  (1962).  Tides,  In  The  Sea,  Vol.  1,  M.N.  Hill,  ed.  John  Wiley  and 
Sons,  New  York,  764-801. 

McLeish,  W. ,  and  D.J.P.  Swift,  R.B.  Long,  D.  Ross  and  G.  Merrill  (1981). 
Ocean  Surface  Patterns  above  Sea-Floor  Bedforms  as  Recorded  by  Radar, 
Southern  Bight  of  North  Sea,  Marine  Geology,  4_3,  MI -MS . 


43 


NOAA  (1979).  Georges  3ank  and  Nantucket  Shoals,  National  Ocean  Survey 
Chart  #13200 . 

NOAA  (1981),  Tidal  Current  Tables  1982,  Atlantic  Coast  of  North  America, 
National  Ocean  Survey  Publ. 

Valenzuela,  G,R.  (1981).  A  Remote  Sensing  Experiment  in  the  Nantucket 
Shoals  (3E3EX),  IUCRM  Symposium  on  ’‘Nave  Dynamics  and  Radio  Probing  of 
the  Ocean  Surface,”  Miami  3each,  Fla.;  submitted  for  pu  lication  in  the 
proceedings,  Plenum  Press,  New  York* 


44 


APPENDIX 


This  appendix  includes  the  computer  programs  developed  for  the  data 
treatment.  The  language  is  BASIC  and  the  documentation  is  self- 
explanatory. 
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PROGRAM  "f-:  I  TEC  fi¬ 
rms  PROGRAM  IS  USED  TO  ENTER 
CURRENT  METER  SPEEDS  AND 
DIRECTIONS  PROM  K  EYBOAPD  AND 
STORE  THEM  ON  TAPE  DATA  MAY  BE 
CORRECTED  AFTER  EVERY  16  ENTRIES 


10  OPTION  BASE  1 
20  DEG 

SO  DIM  Si  16>,  DUS) 

40  C*=“CM2EXP“ 

50  0$=“ CM2DAT" 

60  DISP  “THIS  PROGRAM  STORES  C- 
M  SPEEDS  AND  DIRECTIONS  IN 
FORM  SS.DDD" 

70  DISP  “IN  5  FILES  OF  32  NUMBE 
RS  AN  ADDITIONAL  FILE-’C 
M-EXP'  GIVES” 

SO  DISP  “DATES  START  TIME,  AND 
TIME  INTERVAL” 

SO  CREATE  C * ■ 1 
lOO  CREATE  OT, 5 
110  ASSIGN#  1  TO  D$ 

120  ASSIGN#  2  TO  Ct 
130  C=l 
140  DISP 

150  DISP  "ENTER  DATE  AS  JULIAN  0 
AY" 

160  INPUT  0 

170  DISP  “ENTER  START  TIME, AS  HH 
.  MM” 

1S0  INPUT  T 1 • T2 

190  DISP  “ENTER  TIME  INTERVAL" 
200  INPUT  01 
210  T 1 =Tl *60+T2 
220  PRINT#  2  •  T 1 , D • D 1 
230  FOR  J= 1  TO  10 
240  DISP  “SEGMENT"  ;  J .  “SEQ  NO",C 
250  DISP  “ENTER  16  SPEEDS , D I RECT 
IONS'* 

260  FOR  1=1  TO  16 
270  INPUT  S<I>,D<I> 

2S0  C=C+1 
290  NEXT  I 

300  IMAGE  DO • 5X , ODD , 5X , DUO 
310  DISP 

320  DISP  “NO  ? PD  DIR  “ 

330  FOR  1=1  TO  16 

340  DISP  USING  300  ;  I  •  S <  I  •  0 •  I  .» 
350  NEXT  I 
360  DISP 

370  DISP  "ANY  CORRECTIONS’* 

3S0  INPUT  A* 

390  IF  A S # “ Y E S '*  THEN  430 
400  DISP  “ENTER  LINE  NO  , CORRECT 
SPEEO,  CORRECT  DIRECTION” 
410  INPUT  I  . S (  I  > • D '  I  » 

420  GOTO  320 
430  FOR  1=1  TO  16 
440  D < I )=0< I >'1000 
450  E=S(.  I  >+D<  1  > 

460  PRINT#  1  ;  E 

470  NEXT  I 
480  NEXT  J 

490  DISP  "ENTER  END  DATE  AS  JO" 
500  INPUT  D 

510  DISP  “ENTER  END  TIME, hH . MM" 


5  v  0 

INPUT  T1,T 

530 

T  1  =  T  1*  6  Q  +  T 

tm 

540 

PRINT#  2  ; 

D  >  T 1 

550 

ASSIGN#  1 

TO  4 

560 

ASSIGN#  2 

TO  * 

570 

DISP  "C-M 

OATH 

580 

END 

PROGRAM  "TPNTD  “ 

THIS  PROGRAM  TRANSFERS  CURRENT 
METER  DATA  FROM  TAPE  TO  DISC 


10  OPT  1  UN  W£  ‘ 

20  DEG 

4  0  C $  =  " CM2  EXP • 0 7  0 0 " 

50  D  *  =  "  C  M  2  0  A  T  D 7  0 0 " 

60  DISP  "THIS  PROGRAM  STORE?  - 
M  SPEEDS  AND  DIRECTIONS  IN 
FORM  SS  ODD  ON  01'-''  RROM  TAP 
Ew 

70  DISP  "IN  1  SPILES  OF  32  HL'MPE 
RS  AN  ADDITIONAL  FILE,’C 
M-EXP 1  G I UES " 

ftft  DISP  "ST AFT  TIME-  END  TIME  • 
DATE  HMD  INTER V A L  BETWEEN  D 
AT  A  " 

SO  CREATE  C*,l 
lOO  CREATE  Of >  15 
110  ASSIGN#  1  TO  0* 

120  ASSIGN#  2  TO  CT 
130  T*=”  T “ 

140  DATA  C M4DAT . CM5DAT . CM6DAT 
Ml  ASSIGN#  3  TO  "CM4EXP  T " 

142  READ#  3  ■  A.B 

143  ASSIGN#  3  TO  -‘CMCEXP  1  " 

144  READ#  3  :  F  *  G . C - n  F 

'161  PRINT#  2  *  A  ,  B  *  f.  •  l  n 
162  FOR  J=1  TO  3 

170  READ  FI 

13*  ASSIGN#  3  TO  f $ 

230  FOP  1=1  TO  160 
240  READ#  3  ;  C 
242  PRINT#  1  ,  C 

250  NEXT  I 

2S0  NEXT  J 

320  DISP  "DATA  TPANSFEREDU 

460  END 
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PROGRAM  "VSUM3 


PROGRAM  “  R  i  7  EDRFT 


THIS  PROGRAM  PLOTS  THE 
PROGRESSIVE  ‘-‘ECTOR  DIAGRAM  FOR 
CURRENT  METER  J1  FROM  SPEED  AND 
DIRECTION  DATA  FILED  IN  ‘•CMODAT" 
THROUGH  "CM4DAT" . 

10  OPTION  BASE  1 

11  PLOTTER  IS  705 

20  DEG 

21  XO . Y0=0 

22  DATA  CMODR i , ISO. CM-OAT , 224, C 
M 1  DAT .480  -  CM 2D AT . 430 

23  OATA  CM3DAT. 430,CM4DAT, 144 

31  SCftLE  -30000, 6800 , -260O0 • 100 
0 

32  GCLEAR 

33  MOVE  XO,YO 

34  FUR  J=1  TO  6 

35  REAO  0$,  N 

60  ASSIGN#  1  TO  D$ 

SO  FOP  1=1  TO  N 

90  READ#  1  :  C 

?  1  S  =  I P  <.  C  >  S  D = F  P  <  C  ;•  *  t  O  A  O 
100  X=S*SIN<0>  e  Y  =  S  * C  0 S  »  D .< 

110  X0=X0+X  i$  Y0=Y0+Y 

130  PLOT  X 0 , Y A 

140  PENUP 

150  NEXT  I 

160  NEXT  j 

170  END 


THIS  PROGRAM  CALCULATES  THE  EAST 
AND  NORTH  COMPONENTS  OF  DRIFT  OF 
THE  PROGRESSIVE  VECTOR  DIAGRAM 
OVER  A  25  HOUR  PERIOD  FOR 
CURRENT  METER  X  AND  STORES  THESE 
IN  THE  FORM  OF  SPEEDS  o:m/sec 
IN  FILE  ‘  DR  IFT-X" 

10  OPTION  BASE  1 
20  DEG 

30  DATA  CXI  DAT..  CX2GRT ,  CX3DAT 
40  ASSIGN#  1  TO  "CX0DATM 
50  CREATE  ‘'DRIFT  -v-‘  .  i 
60  ASSIGN#  2  TO  "DRIP1-::" 

70  FOR  J=1  TO  3 
SO  XO , Y0=O 
90  FOR  1=1  TO  200 
1O0  READ#  1  ;  C 

110  3  =  I P 1  C .  '  *5  Q  =  F  P r  C  ■  t  !  0  o  A 
120  X=S*SIN‘0>  6  Y*StlOSvD> 

130  X O  =  X 0 + X  A  Y 0  =  Y Q + Y 
140  IF  J=3  AND  I=80  THEN  GOSUB  2 
20 

150  IF  J  =  5  AND  1  =  160  THEN  GOSUB 
220 

160  IF  J=8  AND  1=40  THEN  GOSUB  2 
20 

170  NEXT  I 

180  X  0 = X  0 ' 2  0  O  {  0  =  Y  0 ' 2  0  0 

190  PRINT#  2  ;  XO,YO 
2O0  NEXT  J 
210  STOP 
220  READ  C * 

230  ASSIGN#  1  TO  Dr 
240  RETURN 
250  END 


47 


PROGRAM  “V3UM8 


IN  THIS  PROGRAM,  THE  PROGRESSIVE 
VECTOR  DIAGRAM  FOR  CURRENT  METER 
X  IS  PLOTTED •  BUT  WITH  THE  DRIFT 
REMOVEO  FOR  EACH  25  HOUR  v  200 
PO IHTS  >  PERIOD  THE  RESULT  IS  ft 
SERIES  OF  OVERLAPPING  TIDAL 
LOOPS  ‘.nidus  ratti>. 

10  OPTION  BASE  1 
20  PLOTTER  IS  705 

30  DEG 

31  SCALE  -1S00, 4000. -2100 ■ 210O 

32  GCLEAP 

40  DATA  CXI  DAT . CX20AT ■ CX30AT 

50  ASSIGN#  1  TO  "CXODAT" 

51  ASSIGN#  2  TO  “DRIFT-X" 

52  PENUP 

5 0  FOR  J = t  TO  9 

51  READ#  2  ,  XI..  VI 
70  X0,Y0=O 

80  FOR  1=1  TO  200 
90  READ#  1  ;  C 
100  S=IP'C>  0  D=FP'C  >T1O00 
120  X»S*SIN'.D->-Xl  ®  Y=S4C03  (.  D  >  -Y 
1 

130  XO=X0+;<  i?  Y0  =  YO+Y 

131  PLOT  X0-Y0 

140  IF  J  =  3  AND  1=80  THEN  GOSUB  3 
40 

150  IF  J=5  AND  1=180  THEN  GOSUB 
340 

180  IF  J=8  AND  1=40  THEN  GOSUB  3 
40 

130  NEXT  I 
210  NEXT  J 
340  READ  D* 

350  ASSIGN#  1  TO  Df  - 
380  RETURN 
370  FRAME 

380  MOVE  -1000.1800 

390  CSIZE  3 •  7.10 

400  LABEL  "CURRENT  METER  X" 

410  CSIZE  2  5 •  7.10 

42a  LABEL  "PROGRESSIVE  VECTOR  01 

|.|  M 

430  LABEL  "OF  SPEED  AND  DIRECT  10 
N“ 

440  LABEL  "FOR  JULY  10  TO  20.193 

450  LABEL  "DRIFT  REMOVED" 

480  MOVE  -1000 . 1200 

470  DRAW  -558, 1200 

471  MOVE  -1000,1000 

480  LABEL  “  2  KM  " 

490  END 
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P  P  0  *  j  P  h  M 


"  V  i;  IJM8  ** 

A  25  HOUR  1  I  DHL  CYCLE  IS 
P  R  0  D  U  C  E  D  6  t  C  A  L  C:  U  L  A  T  l  N  G  T  H  t 
AVERAGE  VAuUES  FOP  POINTS  1  TO 
2 0 0  FOR  t't-t  TEN  Oh Y  CURPEnT 
METER  M  RECORDS  THESE  POINTS 
APE  PLOTTED  AS  PROGRESSIVE 
VECTOR  DIAGRAMS.  THE  25  HOUR 
RVERfl GE3  ARE  THEN  FURTHER 
AVERAGED  INTO  A  12.5  HOUR  <100 
POINT)  TIDAL  CYCLE,  AND  THESE 
POINTS  ARE  STORED  IN  FILE  "MTIDE 

10  OPTION  BASE  1 

11  DIM  ft 2  O  0  >  ,  Y  <  2  0  O  > 

20  PLOTTER  IS  705 

0  DEL 

3 1  SC AL E  - 1 S00 , 4000 , -2 1 O u • 2 1 00 

32  G CLEAR: 

33  FOR  1  =  1  TO  20v? 

34  X  <  1  >  ,  V  \  I  )  -0 

35  NEXT  I 

40  DATA  CM-OAT, CMtDAT, CM2DAT, CM 

3  0  A  T ,  $ 

50  ASSIGN#  t  TO  11 C  MOD  AT " 

51  ASSIGN#  2  TO  "DRIFT-M" 

52  PENUP 

*0  FOR  J=1  TO  9 
El  READ#  2  ;  XI, VI 
SO  FOR  1=1  TO  200 
90  READ#  1  •  C 

10O  S=IP'C  >  G  D»FP<C>*100& 

120  X < I > =  X v I  '  + S * S I H (  D  >  - X 1  3  Y-  l  > 
=  Y K  I  >+S*CQS<D  ' - Y 1 
140  IF  J=1  AND  1  =  1  SO  THEN  GQSUB 
340 

150  IF  J  =  2  AND  1=1 94  THEN  COS US 
340 

ISO  IF  J=5  AND  1=64  THEN  GOSUB  3 

4  0 

170  IF  J=7  AND  1=144  THEN  GOSUB 
340 

130  NEXT  I 

210  NEXT  J 

211  X0,V0=O 

220  FOR  i= 1  TO  200 

230  X<  I  >  =  X<  I  >  '9  (2  V*:  I  >=Yv  I  '  '9 

2  4  A  X  0  =  X  0  + X  I  »  S  V  0  =  Y  0  Y  *  I  > 

250  PLOT  Y.0,Y0 

260  NEXT  1 

270  GOTO  370 

340  READ  DS 

350  ASSIGN#  1  TO  OR 

360  RETURN 

370  FRAME 

330  MOVE  -1O00.1300 

390  CSI2E  3.  7. 10 

400  LABEL  "CURRENT  METER  M" 

410  CSIZE  2  5,  7, 10 
420  LABEL  "PROGRESSIVE  '.'ECTOR  DI 
fiGRRM  OF" 

430  LmBEL  " AVERAGE  SPEED  AND  DIF: 
EOT  ION" 

440  LABEL  "FOP  JULY  10  TO  2019S 

^  II 

450  LABEL  “DRIFT  REMOVED" 

460  M 0 V E  -10O0.  I  2 0 0 

470  DRAW  -556.1200 

471  MOVE  -1060.1003 
430  LABEL  H  2  KM." 

490  CREATE  "MTIDE”,? 

500  ASSIGN#  1  TO  "MTIDE" 

510  FOR  1=1  TO  IDO 


5  2 0  X  =  '•  X  I  >  +  X  *,I  +  1 0  0  >  '  2 
5  3  0  V  =  (.  V I  >  +  Y  <  I  + 1 0  O  *  >  /  Z 
54  0  PRINT#  1  ,  y  .•  V 

550  NEXT  I 
560  END 


P R 0 GRAM  “ 1 h B L E  —  T " 

THIS  PROGRAM  USES  THE  12.5  HOUR 
TIDAL  CYCLE  DATA  STORED  IN 
“MTIDE"  TO  PRINT  A  TIDE  TABLE 
FOR  CURRENT  METER  M 

10  OPTION  BASE  1 
20  DEG 

30  DIM  S 1 2 >  ,  0 '-1 2 > 

40  ASSIGN#  1  TO  "MTIDE" 

50  READ#  1  ;  A, A 

60  PRINT  "  CURRENT  METER  M" 

70  PRINT  "  TIDE  TABLE" 

80  PRINT 

90  PRINT  “TIME  DIRECTION  SPEE 

0" 

10O  PRINT  “hours  decrees  knot 

110  PRINT 

120  FOR  -J=l  TO  12 

130  X0,  r‘0=0 

140  FOR  1=1  TO  8 

15y  READ#  1  ;  X,Y 

160  X0=X0-f X  G  Y0  =  Y0^Y 

170  NEXT  I 

130  X 0  =  X 0 / 3  i?  Y0  =  Y0xS 
190  6=S0P<X0  •2+Y0^2> 

2  O 0  S < J  > =34 . 0194 
210  D < J  >  = A  T  M  2  <  X  0 . Y 0  > 

220  IF  D <  J >  <  0  THEN  D < J> *0 < J 
230  NEXT  J 
240  J=0 

250  PRINT  USING  290  ;  J,.DU2)>^ 
12) 

260  FOR  j=l  TO  11 

270  PRINT  USING  290  ;  J,CKJ>,S<J 
) 

230  NEXT  j 

290  IMAGE  2X , DO , 4X , ODD , 4X  #  0 . DD 
30O  ASSIGN*  1  TO  t 
310  END 
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PRuGRhm  " 7 i DEHOUP" 

THIS  PROGRAM  CALCULATES  AVERAGE 
CURRENT  METER  SPEEDS  HMD 
DIRECTIONS  FOP  EACH  HOUR  'FROM 
t/2  HOUR  BEFORE  TO  l '2  HOUR 
RFTER  THE  HOUR)  AND  LISTS  THESE 
FOR  THE  10  DAYS  OF  OPERATION. 

10  DEG 

20  OPTION  BASE  1 

30  DATA  CN 1  DAT • CN2DAT . CN3BAT 

40  ASSIGN#  1  TO  "CNuDAT" 

41  READ#  1  >  C.C,C,C 

A0  03=  190  a  N=9  a  T=15 

60  PRINT  “TIDE  TABLE  DERIVED  FR 
OM" 

70  PRINT  “CURRENT  METER  ' M '  AT 
DEPTH  OF’.’ 

SO  PRINT  “TWENTY  ONE  METERS" 

90  PRINT  “Speeds  and  directions 
are" 

lOO  PRINT  "averages  ot  eight  mea 
su  remen  t  s “ 

110  PRINT  “  The  table  gives  ho 
url>'  speeds" 

120  PRINT  "and  directions  tor  th 
e  ten  c jaw" 

130  PRINT  •* period  fro m  July  10  t 
o  2  0  .•  1 9  S'  2 " 

140  PRINT  “Hours  are  in  universe 
1  time  ’■ 

L 50  PRINT  “Subtract  tour  hours  i 
or  local" 

tb0  PRINT  “time  or  five  hours  to 
r  Eastern" 

170  PRINT  "Standard  Time  " 

130  PRINT 

130  K  =  1 

2O0  GOSUB  400 

210  FOR  J-i  TO  233 

220  v 0 . Y0=O 

230  FOR  1=1  TO  3 

240  READ#  1  ,  C 

250  3  =  I P <  C  >  *  O=FPo:.)  tcl000 

5b0  X=S*S1N<0>  9  Y=STC0Si:O> 

278  va=X0+X  @  Y0  =  Y0+‘* 

271  IF  J=59  AND  1=4  THEN  GOSUB  5 
10 

272  IF  J=  1 1 9  AND  1=4  THEN  GOSUB 
510 

273  IF  J=*  1  79  ANO  1=4  THEN  GOSUB 
510 

274  IF  J=239  AND  1=4  THEN  540 
280  NEXT  I 

290  X 0  =  X 0 ' 8  9  Y 0  =  Y 0 / 3 

300  S  =  SQR(  X0^2  +  Y0/v2  > 

301  3=S* , 0 l 94 
318  D=flTN2<X0.  Y0') 

320  IF  D<0  THEN  D=D+3b0 
330  PRINT  USING  390  ;  T,S,0 
340  T=T+ 1 

350  IF  T=24  THEN  GOSUB  391 
380  NEXT  J 

390  INAGE  5X,0Q,5X,D  DD , 5X  * ODD 

391  T  =  0 

400  03  =  09"*- 1  9  N  =  N+1 
410  K*K*-1 

420  PRINT  9  PRINT  3  PRINT 


430  IF  K  =  - 1  THEN  PRINT  9  P  R I  N  T  «- 
PRINT  ft  PRINT 


440  PRINT  “ 

“  19821' 

441  PRINT 

453  IF  K=1  THEN 
480  PRINT  “ 

EOT  ION “ 

470  PRINT  41 


DAY  "  09  “  JULY"  ;N 
4  38 

HOUR  SPEED  DIR 
ut  k  t  s .  d e  a 


rees" 

471  PRINT 
480  RETURN 
510  READ  F $ 

520  ASSIGN#  1  TO  F* 
530  RETURN 
540  END 


PROGRAM  "PlICM 


EIGHT  VALUES  OP  SPEED  HMD 
DIRECTION  FOP  CURRENT  METER  X 
ARE  CONVERT  ED  TO  X  BHD  Y  HMD 
AVERAGED,  THEN  RECONVERTED  TO 
GIVE  HOURLY  SPEEDS  AND  DIRECTION 
THESE  ARE  PLOTTED  ACROSS  THREE 
SHEETS  OP  PAPER  PROGRAM  EDITING 
ALLOWS  CURRENT  METERS  M  AND  N  TO 
BE  SIMILARLY  TREATED 


10  DEG 

20  OPTION  BASE  1 

21  PLOTTER  IS  705 

30  DATA  CXI  DAT , CX2DAT, CX3DAT 
40  ASSIGN#  l  TO  “CXODAT" 

50  READ#  1  •;  .C  ,  C  ,  C  .  C 
CO  09=1  SO  ft  N  =  9  ft  T=  1 6 
70  SCALE  -4.100.-100,100 
SO  LINETYPE  1 
90  XAXIS  -100. 24,0, 96 
100  XAXIS  55,24,0,9b 
1 10  K  =0 

120  V  A X  I  S*  0 ,  10,  -  1  OO  .  0 
130  VAX  IS  0,7  5,  lO  ,  100 
140  LINETYPE  3,4 
150  PENUP 

220  FOR  J=1  TO  239 

230  X0, Y0=O 

240  FOR  1=1  TO  0 

250  READ#  1  ;  C 

2b O  3 - I P • C  >  ft  o  =  f P  <  C  > 1 1 0 O 0 

270  X  =  S*SINcD*>  ft  Y=S*COS<D> 

230  XO=X0+X  ft  Y0=YO+Y 
281  IF  J=59  AND  1=4  THEN  GQSUB  5 
70 

290  IF  j  =  1 19  AND  1=4  THEN  GO SUB 
5  70 

300  IF  J  = 1 7 9  AND  1=4  THEN  GQSUB 
570 

330  IF  J  =  239  AND  1=4  THEN  COO 
340  NEXT  I 

350  X 0  =  X  O '8  ft  Y 0  =  Y  0 / 3 
360  S  =  SQP '  X O /N  2  +  Y 0 ^ • 

380  0=A TN2  K  XO  ,  YO 
390  IF  DO  THEN  0=0+368 
400  <  PLOT  J-K.S-100 
410  PLOT  J— K . 0^4+1 0 

420  IF  J=9C  THEN  GOSUB  440 

421  IF  J= 192  THEN  GOSUB  448 
430  NEXT  J 

440  BEEP 

44  1  DISP  '‘CHANGE  PAPER" 

442  DISP  "PRESS  CONTINUE" 

450  PAUSE 

451  LINETYPE  1 

460  XHXIS  -100.24,0-96 

461  XAXIS  55,24,0,96 

462  K  =  K  +  96 

463  VAX  I S  0,10,-100,0 

464  YAXIS  0*7.5,10, 100 

465  LINETYPE  8,4 

466  PENUP 
470  RETURN 
570  READ  Ft 

530  ASSIGN#  1  TO  F$ 

590  RETURN 
CAO  END 


+' P 0 G R AM  * ■  X - N -UP 0 G  J L  " 

THIS  PROGRAM  PLOTS  A  PROGRESS '"F 
VECTOR  DIAGRAM  FOP  CURRENT  ME c 
X  FOR  THE  FIRST  TWENTY  ONE  H'»nf 
FOR  COMPARISON  WITH  THE  TRACK  v* 
DROGUES  DEPLOYED  WEST  OF  PHELPS 
BANK  . 

lO  DEG 

20  PLOTTER  IS  705 
35  ASSIGN#  1  TO  " CXODAT “ 

40  ASSIGN#  2  TO  “DRIFT -X “ 

50  READ#  2  ;  X 1 , Y 1 
CO  READ#  1  ;  C , C , C , C 

7  0  X  0 ,  Y  0  =  0 
80  GCLEhR 

90  SCALE  -600O, 1 4000 * -4500 , 1 050 
0 

too  FOP  1=1  TO  170 

llO  READ#  i  •  C 

120  S  =  I P  * C  *  ft  D  =  FP < C  1  *  t WOO 

1 3 O  X  =  S * S  I  r J  •  0  1  ft  Y  =S  i  C  0 S  •  D  1 

140  ;  :  =  •  X-‘  1  •  ¥ 4  5  ft  Y  = '  v  Y1  •  t 4  5 

1 5  0  X  0  =  X  0  +  X  ft  V  0  -  i  0  + 

ICO  PENUP 
170  PLOT  >,0  ,  Y0 
180  NEXT  I 

190  MO V E  -50O0.-40 0 0 

200  C  ’SIZE  3 ,  7  1 0 

210  LABEL  "CURRENT  METER  X 

220  C  SIZE  2  5  -  7-1 0 
230  LABEL  '‘14  203  Dav  !9l  to  1131 
3  Dav  192" 

240  END 


P  R  0  GRAM  R  iV  i'  •' 

THIS  PROGRAM  F  L  7  :  ^  -c*;- 
P U  r  I  T  I ' J N Fijp  DP»  -'j'JE  **  *  ■  D 

IN  FILE” A-BOY "  >  ON  THE  S  m  M  t 
SChlE  AS  THE  PLOT  FROM  PROGRAM 
“X-W-DROGUE " 

1  0 P T I  0 N  BASE  1 

2 O  SCALE  -5000-1 5 000 .-4060.  11 00 
0 

30  PLOTTER  IS  705 

31  ASSIGN#  1  TO  “A-BOY" 

32  DISP  "WHAT  KIND  LINE  \ 1 -8  ■  " 

33  INPUT  l 

34  LINETYPE  L 

40  FOR  I=i  TO  30 
50  READ#  1  .  X . Y 

5  1  PLOT  ,  V 
60  NEXT  I 
70  BEEP 
80  END 


51 


